Spin gap and magnetic resonance in superconducting BaFei gNio.i As2 
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We use neutron spectroscopy to determine the nature of the magnetic excitations in supercon- 
ducting BaFe1.9Nio.1As2 (Tc = 20 K). Above Tc the excitations are gapless and centered at the 
commensurate antiferromagnetic wave vector of the parent compound, while the intensity exhibits 
a sinusoidal modulation along the c-axis. As the superconducting state is entered a spin gap gradu- 
ally opens, whose magnitude tracks the T-dependence of the superconducting gap observed by angle 
resolved photoemission. Both the spin gap and magnetic resonance energies are temperature and 
wave vector dependent, but their ratio is the same within uncertainties. These results suggest that 
the spin resonance is a singlet-triplet excitation related to electron pairing and superconductivity. 

PACS numbers: 74.25.Ha, 74.70.-b, 78.70.Nx 



The magnetic scattering in the high-transition- 
temperature (high- Tc) copper oxide superconductors is 
characterized by strong spin correlations in the vicinity of 
the antiferromagnetic (AF) wave vector of the magneti- 
cally ordered parent materials, and a spin 'resonant' mag- 
netic excitation whose energy scales with and whose 
intensity develops like the superconducting order param- 
eter [H,!! Js JIl Like the cuprates, the Fe-based supercon- 
ductors a a 3 are derived from electron [1, 0, [l^j Ei or 
hole do ping of their AF long-ranged ordered parent 



compounds ISL 13. IsL TtL 3] and spin fluctuations 
have been postulated as the possible glue for mediating 
the electron pairing for superconductivity [3, 0, 21 1. 
Indeed, the very recent observation of the same type of 
magnetic resonant excitation in the iron-based supercon- 
ductors 22, 23, 24] inexorably links these two high-Tc 
superconductor families together, and strongly suggests 
that the pairing mechanism has a common origin that is 
intimately tied to the magnetic properties. 

An essential step in elucidating the role of magnetism 
in the superconductivity of these materials is then an in- 
depth determination of the energy [E — huj) and wave 
vector {Q) dependence of the low energy ma gne tic scat- 
tering as the superconducting state is formed [H Sim. 
If electrons in the Fe-based superconductors indeed form 
pairs of spin-singlets below Tc as in conventional super- 
conductors [231 and high- Tc copper oxides, there can be 
an energy associated with exciting the spin-singlet into 
the high-energy spin-triplet state, without unbinding the 
electron pairs. In this picture, the Cooper pairs should 
exhibit a wave-vector-independent spin gap with a T de- 
pendence that gradually opens below Tc, much like the 
temperature dependence of the isotropic superconducting 
gap function observed by angle resolved photoemission 



excitation spectrum S{Q,uj). We find that the spin-gap 
does open gradually below Tc, but the gap energy is dis- 
persive rather than wave-vector independent and tracks 
the dispersion of the resonant mode that has been ob- 
served [24] . These results suggest that the resonant mode 
is indeed the spin-singlet to spin-triplet excitation. 



crystals of superconducting 
Tc = 20 K) because these 



spectroscopy (ARPES) experiments [26|, |27|. We have 
used inelastic neutron scattering to probe the wave vec- 
tor and energy dependence of the low energy magnetic 



We chose single 
BaFe1.9Nio.1As2 (with 
samples have excellent superconducting properties 
In the absence of Ni-doping, BaFe2As2 is a nonsuper- 
conducting metal that orders antiferromagnetically with 
a spin structure shown in Fig. la [l^. Because of the 
unit cell doubling along the orthorhombic a-axis and 
c-axis spin arrangement, magnetic Bragg reflections 
occur at wave vectors Q — (1,0,1) and (1,0,3) type 
positions and are absent at Q = (1,0,0) and (1,0,2) 
[l6l , [TtI . [isj . Previous neutron scattering experiments 
on hole-doped Bao.86Ko.4Fe2As2 powder samples [2^ 
and single crystals of BaFei.84Coo.i6As2 (Tc = 22 K) 
[23! have shown that the effect of superconductivity 
is to induce a neutron spin resonance at energies of 
^Sfc^T, remarkably similar to the doping dependence 
of the resonance in hi gh- T c copper oxides 28, [2^ 
and heavy fermions [30l. l3lj. Measurements on single 
crystals of BaFe1.9Nio.1As2 (Tc = 20 K) ^ suggest 
that the resonance actually exhibits dispersion along 
the c-axis, and occurs at distinctively different energies 
at the three-dimensional (3D) AF ordering wave vector 
Q = (1,0,1) and at Q = (1,0,0). We note that in 
the parent materials the spin wave dispersions in the 
Fe-based superconductors are anisotropic and clearly 
3D in nature, as opposed to the purely two-dimensional 
spin wave dispersion on the parent cuprates. For the 
cuprates the spin fluctuations in the superconducting 
regime are again purely 2D [28, 2^, S^l, while the iron- 
based superconductors appear to exhibit anisotropic 3D 
behavior like their parents. 
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FIG. 1: (color online), (a) Schematic diagram of the Fe 
spin structure in the BaFe2As2, which has magnetic Bragg 
peaks a,t Q — (1,0,1), (1,0,3), etc. For our experiment on 
BaFe1.9Nio.1As2, we use the same unit cell for easy compar- 
ison, (b) Temperature dependence of the spin gap as de- 
termined from energy scans (Fig. 3c) and the temperature 
dependence of the scattering at Q = (1,0, 1) (Fig. 3d). The 
solid curve represents the temperature dependence of the BCS 
gap function. (c,d) Schematic of the magnetic response and 
spin gaps at Q = (1,0,0), and (1,0,1), respectively. Mea- 
surements at Q = (1,0,3) showed similar behavior as those 
at Q = (1,0,1). 



FIG. 2: (color online). Examples of constant energy scans 
around the (1, 0, 3) position for E = 1 meV obtained with 
Ef = 5 meV above and below Tc on SPINS, (a) Q-scan along 
the [H, 0, 3] direction for S = 1 meV at 24 K and 2 K. A 
clear peak centered at (1, 0, 3) at 24 K disappears at 2 K, in- 
dicating the opening of a spin gap. (b) Similar scan along the 
[1, 0, H] direction showing a peak centered at (1, 0, 3) that dis- 
appears below Tc. (c) Using scattering at 2 K as background 
scattering, we determine the normal state L-modulation of 
the spin fluctuations by subtracting the 2 K data from 24 K 
data. It is clear that spin fluctuations are 3D and have sim- 
ilar modulations along the c-axis as spin waves, d) Q-scan 
in the superconducting state through the magnetic resonance 
position, and above Tc near (1,0, 1). 



The neutron scattering measurements were carried 
out on the SPINS cold and BT-7 thermal triple-axis 
spectrometers at the NIST Center for Neutron Re- 
search. We label the momentum transfer Q —{qx, <ly,<lz) 
as {H,K,L) = {qxa/2TT,qyb/2TT,qzc/2TT) reciprocal lat- 
tice units (rlu) using the orthorhombic magnetic unit cell 
of the parent undoped compound (space group Fnimm, 
a = 5.564, h — 5.564, and c — 12.77 A) for easy com- 
parison with previous spin wave measurements on the 
parent compounds, even though the actual crystal struc- 
ture is tetragonal 33, 34,[35|. Many single crystals were 
co-aligned to obtain a total mass of ~1.2 grams. The 
in-plane and out-of-plane mosaics of the aligned crys- 
tal assembly are 1.3° and 4.3° full width at half max- 



imum (FWHM), respectively [2J]. For the experiment, 
the BaFe1.9Nio.1As2 crystal assembly was mounted in the 
[H, 0, L] zone inside a liquid He cryostat. The final neu- 
tron wave vector was fixed at either Ef = 5 meV with a 
cold Be filter or at Ef = 14.7 meV with a PG filter in 
front of the analyzer. 

We first probe the wave vector dependence of the low- 
energy spin fluctuations. Figures 2a and 2b show [H, 0, 3] 
and [1, 0, L] scans at -B = 1 meV through the 3D (1, 0, 3) 
Bragg peak position below and above Tc. We see that 
the spin excitations observed above Tc vanish at low T. 



Fig. 2c shows the intensity of the scattering above Tc as 
a function of wave vector along the c-axis, using the low 
T data as background, and reveals the intrinsic wave vec- 
tor modulation of the intensity of the normal state spin 
fluctuations. The solid curve is a flt to the data using 
AS'(Q,cj)(24 K-2 K)= AF{Q)^ sin^{TTL/2) + C, where 
F{Q) is the magnetic form factor of Fe^+ and C is con- 
stant. These data are consistent with previous work on 
BaFe1.9Nio.1As2 which showed that the spin fluctuation 
intensity has a c-axis modulation at E = 8.5 meV, and 
a gap in the superconducting state 2^. For compari- 
son, Fig. 2d shows the magnetic scattering through the 
[1,0, 1] position in the superconducting state at the reso- 
nance energy oi E = 7 meV, and the magnetic scattering 
above Tc- We note that in the undoped AF state, the 
spin wave spectrum in BaFe2As2 has a gap of 9.8 meV 
35], while in the normal state of the doped system we 
find that the spin fluctuation spectrum is gapless. 

The behavior of the low energy spin excitations as a 
function of temperature is shown in Fig. 3, which sum- 
marizes the BT-7 and SPINS data around Q = (1,0, 1). 
Figure 3a shows wave vector [H, 0, 1] scans through the 
Q — (1,0,1) position above and below Tc at E = 2 
meV. A clear Gaussian peak centered at Q = (1,0,1) in 
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FIG. 3: (color online), (a) Q-scans at the [iif, 0, 1] direction 
above and below Tc at hco = 2 meV. The data show the open- 
ing of a spin gap at 2 meV below Tc- (b) Q-scans along the 
[1,0, L] (signal) and [1.2,0,1/] (background) positions, show- 
ing the L modulation of the intensity, with maxima at (1, 0, 1) 
and (1, 0, 3). (c) Constant-(3 scans at Q = (1, 0, 1) at various 
temperatures. The differences between low and high temper- 
ature data show negative scattering due to the opening of a 
spin gap. The data suggest a spin gap value of 1.5 meV at 15 
K and 3.0 meV at 2 K. (d) Temperature dependence of the 
scattering at Q — (1,0, 1) and E — 1 meV shows a sudden 
drop below 18 K (= Tc - 2 K) suggesting that the £ = 1 meV 
spin gap opens at a temperature slightly below Tc. 



the normal state vanishes below Tc, demonstrating that 
the spin gap Agpin > 2 meV. Figure 3b plots the signal 
and background scattering along the [1, 0, L] direction for 
E = 2 meV at 30 K, where we find that the normal state 
scattering also peaks at 3D AF wave vector positions. To 
determine the spin gap value at Q = (1,0, 1), we carried 
out temperature dependent measurements at 2 K, 15 K, 
and 30 K using SPINS. We find a clear reduction in scat- 
tering (net negative values in the subtraction) below 3 
meV and 1.5 meV at 2 K and 15 K, respectively. These 
results show that the maximum magnitude of the spin 
gap at the Q = (1,0,1) wave vector is 3 meV, and the 
energy gap is temperature dependent. 

To quantitatively determine the wave vector depen- 
dence of the spin gap in the superconducting state, we 
carried out constant-Q scans at the Q = (1,0,0), and 
(1,0,1) wave vectors, and collected background data at 
Q = (1.2,0,0), (1.2,0,1), above and below T^ (Figs. 4a 
and 4b). In the normal state (open circles) the mag- 
netic scattering above background at both Q ~ (1,0,0) 
and (1,0,1) appears to increase with decreasing energy 
near the elastic line, and thus suggests that this com- 
ponent of the scattering is quasielastic in nature (peaks 
at E = 0). In the superconducting state (solid squares), 
the low energy scattering is suppressed, while the higher- 
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FIG. 4: (color online). Constant-(5 scans around the (a) 
Q = (1,0,0) and (b) (1,0,1) positions above and below Tc, 
showing the developement of the spin gap at low energies, 
and the enhancement of the magnetic scattering at the res- 
onance energy at each wave vector. Background data are 
indicated by the open diamonds and dashed curves, and were 
collected at Q = (1.2,0,0) and Q = (1.2,0,1), respectively, 
(c, d) x"{Qj^) above and below Tc obtained by subtract- 
ing background and removing the thermal factor (see text). 
Also shown are values obtained from the Q-scans at various 
energies above and below Tc. At both wave vectors there 
is a clear magnetic intensity gain at the resonace energy of 
E = 9.0 meV at g = (1, 0, 0) and 7 meV at Q = (1, 0, 1), and 
a spin gap of 4.3 and 2.5 meV, respectively. The solid lines 
show fits using the model described in the text. 



energy scattering increases in intensity. The overall be- 
havior of the data is remarkably similar to that in the 
optimally hole-doped Lai .T;Sr3;Cu04 [32| and electron- 
doped Ndi.85Ceo.i5Cu04 [29|. However, it is also clear 
that the spin-gap occurs at a lower energy at Q = (1, 0, 1) 
than for Q = (1,0,0), which is quite different than the 
cuprates [34]. Figure 4c, d presents the data in the form 
of the dynamic susceptibility ^"(Q,^), which is related 
to S(Q, ut) through the (removal of the) detailed balance 
factor; x"(Q,w) = {1 - exp{-huj/kBT))S{Q,uj). Re- 
call that the thermal population factor increases with 
decreasing temperature, and this function is divided into 
S{Q,uj) to obtain x"(Q,w) [with x"(Q, w = 0) = 0]. The 
filled circles are x"((5, to) obtained from Q-scans as a con- 
sistency check. Upon entering the superconducting state, 
the spectral weight is rearranged, with the suppression of 
low energy spin fluctuations and the appearance of the 
neutron spin resonance at energies above the spin gap. 
The present data are consistent with the reported spin 
resonance values oi E = 9 meV for Q = (1,0,0) and 
E = 7 meV for Q = (1,0, 1) 24]. We estimate that the 
intensity of the resonance is approximately compensated 
by the opening of the spin gap below the resonance. 
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To quantify the magnitude of the spin gaps at Q = 
(1,0, 0) and (1, 0, 1) in the superconducting state, we fol- 
low previous work [s^l and fit the data with 



S{Q,u;) 



AET 



(r2 + (?ia;)2)(l - exp(-?iw/fcBr)) ' 



(1) 



where E' = Re[{huj - + iV s){nuj + + iTs)Y^'^, A is 
the amplitude, A is the spin gap, T is the inverse life- 
time of the spin fluctuations with fiuj 3> A, i?' is an odd 
function of ii' = hu, and Ts is the inverse lifetime of the 
fluctuations at the gap edge. The solid curves are the 
results of these fits. In the normal state, this functional 
form docs not provide an adequate fit over the entire 
energy range, and we restricted it to lower energies (as 
indicated by the extent of the curve for those data). We 
find A = for both Q = (1,0,0) and (1,0, 1). On cooling 
into the superconducting state, Eq. (1) can be used over 
the entire energy range of the data, and the least-squares 
fit to the Q = (1, 0, 0) data (solid curves in Figs. 4a and 
4c) yields A = 56.7±7.9, T = 13±6.5 meV, A = 4.3±0.8 
meV, = ± 0.73 meV. Similarly, for Q = (1,0, 1) we 
find A = 55.5 ± 14.5, T = 5 ± 0.7 meV, A = 2.5 ± 0.08 
meV, Fs = ± 0.53 meV (solid curves in Figs. 4b and 
4d) . The results of this analysis show that the the super- 
conducting spin gap values for Q = (1,0,0) and (1,0, 1) 
are distinctively different. 

The present measurements, as well as the previous data 
on this material 2J] , demonstrate that the resonance oc- 
curs at i? = 9 meV for Q = (1,0,0), which has a spin 
gap A = 4.3 ± 0.8 meV. For Q = (1, 0, 1) the resonance 
is at the lower energy oi E — 7 meV, and the spin gap 
also occurs at the lower energy of A = 2.5 ± 0.08 meV. 
Therefore these two energy scales track one another, with 
a ratio that is the same within the uncertainties of the ex- 
periments. This is the expected behavior for the singlet- 
triplet transition of a Cooper pair [36| . 

We summarize in Figs. Ib-ld the key results of our ex- 
periments. The measured temperature dependence of the 
spin gap at Q = (1,0, 1) is shown in Fig. lb. The solid 
curve shows the prediction of a simple BCS gap function 
near T^, A(r) = A{1 - (T/Tc)f/'^, which describes the 
data fairly well. Figures Ic and Id plot schematically 
the spin gap and resonance at Q = (1, 0, 0) and (1, 0, 1). 
The two energies exhibit the same dependence on wave 
vector. In ARPES experiments 27|, two isotropic su- 
perconducting gaps with values of 7 meV and 4.5 meV 
were observed for BaFei.85Coo.i5As2 with = 25.5 K. 
Comparison with the Q = (1, 0, 0) neutron measurements 
suggests that the resonance energy at Q = (1,0,0) is 
indeed less than twice the superconducting gap energy. 
These results are consistent with the idea that the reso- 



nance is a bond state related to singlet-triplet excitations 
of Cooper pairs, with a superconducting gap that varies 
with the momentum transfer along the c-axis (scj . 

We thank Songxue Chi, Jun Zhao, and Leland Har- 
riger for coaligning some of the single crystals used in 
the present experiment. This work is supported by the 
U.S. DOE BES No. DE-FG02-05ER46202, NSF DMR- 
0756568, and in part by the U.S. DOE, Division of Sci- 
entific User Facilities. The work at Zhejiang University 
is supported by the NSF of China. 



[1 
[2: 

[s: 
[4: 
[s: 
[6: 
[7: 
[s: 
[9: 

[10: 
[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 

[21 
[22; 

[23' 
[24' 
[25| 

[26' 
[27 
[28' 
[29 
[30^ 
[31 
[32' 
[33' 
[34; 

[35' 
[36 



Electronic address: daip@ornl.gov| 

J. Rossat-Mignod et ai, Physica C 185, 86 (1991). 

H. A. Mook et al, Phys. Rev. Lett. 70, 3490 (1993). 
Pengcheng Dai et al, Phys. Rev. B 63, 054525 (2001). 
C. Stock et ai, Phys. Rev. B 69, 014502 (2005). 

C. K. Xu and S. Sachdev, Nat. Phys. 4, 898 (2008). 

M. R. Norman, Physics 1, 21 (2008). 

S. A. Kivelson and H. Yao, Nat. Mater. 7, 927 (2008). 

Y. Kamihara et al, J. Am. Chem. Soc. 130, 3296 (2008). 

A. Leithe- Jasper et ai, Phys. Rev. Lett. 101, 207004 

(2008). 

A. S. Sefat et al, Phys. Rev. Lett. 101, 117004 (2008). 

L. J. Li et al, arXiv: 0809.2009vl. 

M. Rotter et ai, Phys. Rev. Lett. 101, 107006 (2008). 

C. de la Cruz et ai, Nature (London) 453, 899 (2008). 

M. A. McGuire et al, Phys. Rev. B 78, 094517 (2008). 

J. Zhao et ai, Nature Materials 7, 953 (2008). 

Q. Huang et al, Phys. Rev. Lett. 101, 257003 (2008). 

J. Zhao et al, Phys. Rev. B 78, 140504(R) (2008). 

A. I. Goldman et ai, Phys. Rev. B 78, 100506(R) (2008). 

I. I. Mazin and J. Schmalian, arXiv:0901.4790vl. 

A. V. Chubukov, D. V. Efremov, and I. Eremin, Phys. 
Rev. B 78, 134512 (2008). 

V. Stanev et al, Phys. Rev. B 78, 184509 (2008). 

A. D. Christianson et al, Nature (London) 456, 930 

(2008). 

M. D. Lumsden et al, arXiv:0811.4755 1'l. 

Songxue Chi et al, arXiv:0812.1354vl. 

J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. 

Rev. 61, 1175 (1957). 

H. Ding et al, Europ hys. Lett. 83, 47001 (2008). 
Terashima et al, arXiv:0812.3704vl. 
D. Wilson et al, Nature (London) 442, 59 (2006). 
Zhao et al, Phys. Rev. Lett. 99, 017001 (2007). 
Metoki et al, Phys. Rev. Lett. 80, 5417 (1998). 
Stock et al, Phys. Rev. Lett. 100, 087001 (2008). 
Lake et al. Nature (London) 400, 43 (1999). 
Zhao et al, Phys. Rev. Lett. 101, 167203 (2008). 
J. McQueeney et al, Phys. Rev. Lett. 101, 227205 
(2008). 

K. Matan et al, arXiv: 0810.4790vl. 
M. Eschrig, Adv. Phys. 55, 47 (2006). 



